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We obtain a minimal form of the two-derivative three-nucleon contact Lagrangian, by imposing 
all constraints deriving from discrete symmetries, Fierz identities and Poincare covariance. The 
resulting interaction, depending on 10 unknown low-energy constants, leads to a three-nucleon 
potential which we give in a local form in configuration space. We also consider the leading (no- 
derivative) four-nucleon interaction and show that there exists only one independent operator. 

T— I ■ 
T— I ■ 

O ! I. INTRODUCTION 

(N . 

Ch ■ Purely contact interactions are crucial components of the nucleon-nucleon (NN) and multi-nucleon forces as derived 
^ ^ ] in chiral effective field theory (xEFT) [H-lEl • They encode the short distance properties of the nuclear interaction as 
opposed to terms involving pion exchanges, which are of larger range. At sufficiently low energy, even the pion can 

■ be integrated out, giving rise to the pionless effective theory. The contact vertices are the same in both versions of 
, the theory, the only difference being in the value of the accompanying low-energy constants: in the pionless theory 

I— i' the latter implicitly include the effect of pions, considered as heavy particles. At the level of the NN system two 
, purely contact terms appear at leading order (LO) and contribute to the central part of the interaction p-[ll|- In the 
I ■ three-nucleon (3iV) s yste m the first non-vanishing contribution appears at next-to-next-to-leading order (N2L0) in 
' Q ' the chiral expansion [ill . Il2j . This term has already been included in some three-nucleon interactions (TNI) as the 
^ . i?-term [l^ [l3|. Indeed the inclusion of such term was found to be mandated, in the framework of pionless EFT, 
^ ' by the requirement of renormalizability [T5l - [T8| . In the present paper we focus on the subleading 37V contact terms, 
i.e. those containing two space-derivatives of nucleon fields, and stick to isospin symmetric operators. Only one such 
' term was found to be necessary in Ref. [l^ in order to achieve cutoff independence (see however Refs. [lo, [2l| for a 
^ ' different claim). We would like to stress however, that in the EFT framework one has to include all terms allowed 

■ by symmetry at a given order, not just the one needed by the renormalization procedure, since the EFT is the most 
0^ ' general theory encoding the given symmetry properties. The plan of the paper is the following: in Section |ll] we 

illustrate our strategy to determine a minimal set of 3iV contact operators by considering the leading non-derivative 
operators, thus reobtaining the result that only one operator arises at this order; in Section IlIII the same strategy is 
. applied to the list of all possible two-derivative operators, thus reducing their number from 146 down to 14; further 

■ constraints from relativity, as pertinent to momentum-dependent interactions, are discussed in Section [IVI and allow 
] . to reduce the number of independent operators to 10; we give in Section |V] the resulting 3iV potential in coordinate 

. . ' space in a local form, by choosing an appropriate momentum cutoff. Finally, in Section IVII we consider the leading 
^ ' four- nucleon (4iV) contact Lagrangian and find that there exists only one such operator. 

j3 , II. FIERZ CONSTRAINTS ON THE LEADING CONTACT INTERACTION 

Rotational, isospin and time-reversal invariance allows to list 6 possible operators describing three-nucleon contact 
interactions without derivatives. Indeed, since A^^crTV is odd under time-reversal, an even (odd) number of a matrices 
has to be associated with a p urely real (imaginary) isospin structure. Therefore, the leading order 3iV contact 
Lagrangian density reads |12| : 

= -EiN''NN''NN''N - E2N''a'NN''a'NN''N 

-EsN^a^NN^aWNN^T^N ~ Ef,e'^^ e''^" N'^ N N"^ N t" N 

- -Ei^.of (1) 
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Symmetry properties under permutations are encoded in the anticommuting nature of the nucleon field N. The 
properties under exchange of spin indeces can be conveniently expressed as 

= + i(l][c7^) - ^e^^''ia^][<j'^) (2) 

where 1 is the identity operator in the one-particle spin space, and (,) and [,] denote spin indeces of the enclosed 
operator. Similar relations hold in the one-particle isospin space for the identity and r* operators. Simultaneous 
Fierz rearrangements of spin and isospin indeces of nucleon fields 1-2, 1-3 and 2-3 allow then to derive a set of linear 
relations among the above operators, in a similar way as was done for the (parity violating) two-nucleon contact 
interaction in Ref. [23| . For instance, from the permutation of nucleons 2-3 we find, 
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leaving only one independent operator (in agreement with Ref. |12|). The relations deriving from the remaining 

permutations are linearly dependent on the ones above. The remaining operator was chosen to be of in Ref. [l 
and it gives rise, in momentum space to a contact potential 



with E denoting the corresponding low-energy constant (LEG). Choosing a cutoff depending on momentum transfer, 
it is possible to obtain the corresponding coordinate space potential in a local form [13]. A sensitivity study of bound 
state and scattering observables in ^ < 4 systems to this and other components of the TNI was performed in Ref. [22| ■ 
It should be noticed that the equivalence in the choice of the contact operator is true at the Lagrangian level, but it 
is in general spoiled at the level of the 3iV potential by the cut-off, which may involve non-symmetrical combinations 
of the nucleon momenta. At any rate, in the effective theory framework, different choices are equivalent up to higher 
order corrections. 



III. FIERZ CONSTRAINTS ON THE SUBLEADING CONTACT INTERACTION 



Parity requires that the next-to-leading order 3A'^ contact Lagrangian contain two spatial derivatives. Using trans- 
lational invariance (or momentum conservation) the only possible space-structures can be taken to be of the form 

^L,- V. iV)(iVt V, 7V)(7Vt7V) 

= ^» N){N^N){N^N), 

where the i is required by the hermiticity condition and (hereafter) the superscripts denote the time-reversal properties. 
The relevant isospin structures are 

= 1, Tl • T2, Tl ■ T3, T2 ■ X3 , Xl X T2 • T3, (6) 

where the subscripts of the Pauli matrices refer to the nucleon bilinears they belong to. Even (odd) combinations 
of X (g) T structures under time-reversal have to be associated with spin structures containing even (odd) numbers 
of <T matrices. Finally, the spin-space indeces have to be contracted with Kronecker 5's or Levi-Civita tensors e's. 
Following the above procedure we have obtained a list of 146 operators, displayed in Table HI Fierz rearrangements 
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V^l ■ V2[l, T-i ■ -Ta, T-i ■ 1-3] 




073 




i Vl ■ 1^3 V2 • 0^2(1-1 


X T2 • 


T3I 






04- 


6 


V"l • 0*'l '*^2 ■ [1, Tl • •'"2, Xi ■ -7-3] 




074 




i Vl ■ V2 0*2 ■ ~^3 [1-1 


X T2 • 


T3l 






07- 


9 


V'l • '^■2 ^ 1 ■ ~^ 1 [1, T\ ■ T-l.Ti ■ T3] 




075- 


78 


i Vl X V2 ■ 0^1 [1, 1-1 • 1-2, Tl ■ 


T3, 


T2 • 


T3] 




010 


-12 


Vl • "^2 0*'l ■ ~(^2 [1, Tl • 1-2, -Ti ■ 1-3] 




079- 


82 


i Vl X V2 ■ 0*'2[1, 1-1 • 1-2, 


Tl ■ 


T3, 


T2 ■ 


T3] 




013 


-16 


V'l • 1^1 V2 ■ 0^3[1, Tl ■ ''"2, Ti • 1-3, 1-2 ■ 




083- 


86 


i Vl X V2 ■ 1^3 [1, 1-1 • 1-2, 


Tl ■ 


T3, 


T2 ■ 


T3] 




On 


-20 


Vl • 0^3 ^2 ■ 0^l[l, Tl ■ T-2. 1-1 • 1-3, 1-2 ■ 


■rs] 


087- 


90 


i Vl ■ V2 ~<?1 X 1^2 ■ 0^3(1, Tl 


T2 


Tl 


■ T3 


T2 


T3] 


Oil 


-24 


Vi • ^2 ■ "^3[1, 1-1 ■ T2, 1-1 -1-3,1-2 ■ 


^3] 


O9I- 


94 


i Vl ■ "^1 V2 X ■S'2 • 0^3[1> Tl 


T2,Ti 


■ T3 


T2 


T3] 


025 




V'l X V2 ■ ~^i[ti X 1-2 ■ 1-3] 




095- 


98 


i Vl ■ 0^2 V2 X ■ff'i • ■S'3[l, Tl 


T2 


Tl 


■ T3 


T2 


T3] 


026 




Vl X V2 ■ 0^3 [ti X 1-2 ■ 1-3] 




099- 


102 


i Vl ■ 1^3 V2 X ■ "(7'2[1, Tl 


T2 


Tl 


■ T3 


T2 


T3] 


027 




Vl • V2 O^l X 'a'2 ■ '(^al'Tl X 1-2 ■ X3 




0103 


-106 


j Vl X 1?'2 ■ '5*'3 V2 • O^l [1, Tl 


T2 


Tl 


■ T3 


T2 


T3] 


028 




Vl • O^i 0^2 X 0^3 ■ V2[l-1 X T2 ■ 1-3 




O107 


-110 


i Vl X O^i ■ "^■3 V2 ■ ~<^2l'i-; Tl 


T2 


Tl 


■ T3 


T2 


T3] 


029 




Vl • ■?'2 0^1 X 0^3 ■ Vafxi X T2 ■ 1-3 




Olll 


-114 


i Vl X 1^1 ■ 1^2 V2 ■ 0*'3[1, Tl 


T2 


Tl 


■ T3 


T2 


T3] 


O30 




Vl • 0^3 0^1 X 0*'2 ■ Vafl-l X 1-2 ■ T3 




0115 


-118 


i Vl X V2 ■ "S'l 0^2 • 0^3 [1. Tl 


T2 


Tl 


■ T3 , T2 


T3] 


031 




Vl X V2 • "^1 0*'2 ■ "^3(1-1 X 1-2 ■ X3 




0119 


-122 


i Vl X V2 ■ '^2 0^1 • 0^3[1, Tl 


T2 


Tl 


■ T3 


T2 


T3] 


032 




Vl X V2 • 0*'3 0*'l ■ "^2(1-1 X 1-2 ■ X3 




0123 


-126 


i Vl X V2 ■ '^3 0^1 • 0^2(1, Tl 


T2 


Tl 


■ T3 , T2 


T3] 


033 


-64 


same as 01-32 with V — > V 




0127 


-129 


Vl ■ Vl[l,Ti • T2 


T2 


T3 








065 




i Vl ■ VafTi X T2 ■ T3] 




0130 


-133 


Vl ■ Vl I?! ■ 1?2 [1, Tl ■ T2 


, T2 


■ T3 


I Tl 


•T3 




066 




i Vl • "0^1 V2 ■ ~<^2 [1-1 X X2 ■ X3] 




O134 


-136 


Vl ■ Vl ■ff'2 ■ a^3[l,Tl 


■ T2 


, T2 


■ T3 






067 




i Vl • 0*'2 V2 ■ 0*'i [ti X X2 ■ X3] 




O137 


-140 


V 1 I V 1 • 1?'2 [1 , Tl • T2 , T2 


■ T3 


, Tl 


•T3 




068 




i Vl • V2 0^1 ■ 0^2 [1-1 X T2 ■ 1-3] 




0141 


-143 


Vl - -^2 Vl ■^3[l,Tl 


■ T2 


, T2 


■ T3 


] 




069 




i Vl • 1^1 V2 ■ 0^3 [1-1 X 1-2 ■ X3] 




0144 




Vl ■ Vl ~^lX~^2- ■^3[1- 


Tl 


T2 


T2 


T3] 




070 




i Vl • 1^3 V2 ■ 0^1 [ti X 1-2 ■ T3] 




0145 




Vl ■ 0*1 Vl X "^2 • ■^3[1- 


Tl 


T2 


T2 


T3] 




071 




i Vl • V2 0^1 ■ 0^3 [1-1 X 1-2 ■ T3] 




0146 




Vl ■ 0*2 Vl X "^1 • ■^3[1- 


Tl 


T2 


T2 


T3] 




072 




i Vl • 1^2 V2 ■ 1^3 [1-1 X 1-2 ■ X3] 





















TABLE I: Complete list of 2-derivative three-nucleon contact operators compatible with rotational, isospin, parity and 
time-reversal invariance. Subscripts refer to the nucleon bilinear on which the operators act. For instance 073 = 



also involve the field derivatives: for instance, under permutation of nucleons 1-2, 

Vl ^ i(^2 + V2-'^l+ Vl), "^1 i(^2 + Va+'^i- Vl). (7) 

Out of the 3 X 146 relations from the Fierz identities, 132 are linearly independent, and displayed in the Appendix El 
There remain 14 independent operators, which can be chosen to be 

Oi, 02, 04, 05, O33, O34, O35, 036, O37, 040, 042, 043, O45, 058. (8) 



IV. FURTHER CONSTRAINTS FROM RELATIVITY 



We still have to impose the requirements of Poincare covariance. They can be implemented order by order in the 
low-energy expansion, as detailed in Ref. [24]. As a result, the subleading 37V effective Hamiltonian consists of a set 
of terms whose strength is fixed by the lowest order coupling constant, and a set of free terms which have to commute 
with the lowest order free boost operator. In the following we focus our attention on the latter set. With the choice 
for the non-relativistic nucleon field N{x), 



N{x) 



dp 



bs{p) Xs e" 



(9) 



where bs{p) and 6j(p) are annihilation and creation operators for a nucleon in spin state s, satisfying standard 

bs{p) , bl, (p') = (27r)^(5(p — p') Sss' , the leading order free boost operators Kq acts 



anticommutation relations, i.e. 
in the following way [l^l 

[Ko , 63 (p)] = -imVp bs{p) . 

Only the first 4, out of the 14 operators listed in Eq. ([8]), do not commute with Kq 

[Ko,Oi] = -4m Vl, [Ko,02] ^ 

[Ko,04] -4m Vl • ■(T'l 0^2, [Ko,05] 



4m V 1 Tl ■ T2, 

4m Vi.^1^2 



(10) 



(11) 
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in the same notations as in Table ID The vector operators in the right hand sides of Eq. pT|) . in turn, mix with 45 
others under Fierz rearrangements; for instance, by exchanging the indeces of particles 1-2 we get the identity 

^1 • 1^1 — T ^1 +T ^1 T\- T2 7 V"l 0^1 • 0^2 -T ^1 ~^ I ' ~^ 2 T\ ■ T2 

4 4 4 4 

+ - Vi ■ 0^1 7^2 +^ V"i • 0^1 0^2 Tl • T2 + - Vi • a^2 O^l +^ ■ 0^2 0^1 Ti ■ X2 

X 0^2-^j'^l X a^2Ti •X2 - l*"^! X 0^2 T2 • T3 . (12) 

By analyzing all the 3x49 relations we find that the four operators on the right hand sides of Eq. (jlip are linearly 
independent and, in fact, form a basis of all vector operators satisfying the symmetry requirements. Therefore 
we conclude that the operators 01^2,4,5 are forbidden by Poincare symmetry, and we are left with the remaining 10 
operators. The same conclusion is reached by starting the analysis with an operator basis written in terms of gradients 
with respect to relative coordinates (or Jacobi momenta, in conjugate space), 

VaOcVa-Vi, Vfecx V3-^(Vi + V2), (13) 

in which case one can write an initial list of 116 operators and reduce their number to 10 using the Fierz's constraints. 
This in turns provides a non-trivial check of our calculation. A minimal form of the 2-derivative 3 A'^ contact Lagrangian 
is given by 

= E['^{N'^N)-^{N'^N)N'^N + E'2'^{N'^t''N)-'^{NK''N)N'^N 

+E's ^{N^'t'^N) ■ ^{N''N)N''t''N + E'^^ ■{N''-^ N) ^■{N''^ N)N''N 

+E'^ -(iVt t'^N) ^ ■{N^~^ t''N)N^N + E'e i{N^ N) j{N^a'T''N)N'< N 

+E'^ "^.(iVV^'iV) ^,{N'^a^N)N'^N + E'^ ,{N'' t" N) '^,(iVVViV)iVl'iV 

+E'q ^ -(iVt N) ^{N^N) ■N^-^N + E'^qe"^" ^{N^t^N) x ^(TVV^A) • N'' t^N. 

(14) 



V. THE SUBLEADING THREE-NUCLEON CONTACT POTENTIAL 

The 3A'' potential is obtained by taking the matrix element of the interaction between 3A'' states. Denoting 
ki = Pi — p'^ and Qi = Pi -I- p^, Pi and p,^ being the initial and final momenta of nucleon i, the potential in momentum 
space is found to be 



with 



-Eikf - i?2kj Ti ■ Tj - Eskfcr, ■ a-j - Eikfcr, ■ a-jn ■ Tj 

-i?5 (3ki • (Tiki • (Tj - kf) - Eg (3ki • cr^ki ■ crj - ■ Tj 

--Erh, X (Qj - Qj) ■ (cr, + a-j) + -Egk, x (Q, - Q^) • (cr, + (Tj)Tj ■ 



--E-gki • cTikj ■ (Tj — -Bioki ■ cTikj • cr 



(15) 



El = 


\E[ 


1 3 771/ 3 T-i/ 

+ 2^2 - 2^3 


- \E'j 


+ 1^8' 


E2 = 


E2, 








E3 = 


Ey 


- S3 + \E'q, 






Ei^ 


hE'2 


1 771/ 1 771/ 
- 3^3 - 3^7 


+ E's, 




E5 = 


E'2- 


- E'^- E'-j + \E'^, 




Ee = 


^E'2 


1 171/ 1 ZPl 
^ 3^3 - 3^7' 






Et = 


6E^ 


- QE'^ - 12E'q 


- &E'j 


~^E[, 


Es = 


2E', 


- 2i?^ - 4i?^ - 


- 2E'j - 




Eq^ 




— 3-^3 — -E4 — 


3E't + 


Eg, 


Eio = 


-E'2 


- E!j - E'^ - E'q - E!, 





(16) 
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A local form of the 3N potential in configuration space can be obtained by using a momentum cutoff depending 
only on k^, as done in Ref. [l^, e.g. P{^^: A)F(k^.; A). In this way the result of the Fourier transform is expressed in 
terms of the function 

^o(r;A)= I ^e'P '-i^(p2;A) (17) 
and derivatives theoreof. Explicitly, omitting the argument A in the function Zq, 

F = ^ {Ei+ E2Ti ■ Tj + E'sCTj ■ (Tj + E^Ti ■ Tj(Ti ■ (Tj) 

+{E5 + EeT,-Tj)Sij 
+{E7 + E^T, ■ Xfc)(L • 5),,^^Zo(r,fe) 

+ {Eq + EiqTj ■ Tk)(Tj ■ TijCTk ■ rikZQ{rij)Zlj{rik) (18) 

where Sij and (L • S)ij are respectively the tensor and spin-orbit operators for particles i and j. Notice that a choice 
of basis has been made such that most terms in the potential can be viewed as an ordinary interaction of particles 
ij with a further dependence on the coordinate of the third particle. In particular the terms multiplying Ey and Eg 
are of a spin-orbit character. A specific combination of both has been suggested in Ref. [25*1 . We also observe that 
some of the spin-isospin structures implied by Eq. ([T}, which were equivalent up to cutoff effects, are resolved at the 
two-derivative level. 



Zoink) 



Za{rij) 



Zoirij) 



Zoink) 



VI. LEADING FOUR-NUCLEON LAGRANGIAN 

From the previous discussion, the 3A^ contact interaction consists of a leading contribution, at the N2L0 of the 
chiral expansion, and a subleading one, arising at the N4L0. Parity requires that no SN contact interactions appear 
at N3L0 or N5L0. On the other hand, the leading 4iV contact interaction arises at N5L0, therefore one is lead to 
consider such terms at the same level of accuracy, at least in the framework of pionless EFT. By listing all possible 
spin-isospin structures which respect the discrete symmetries of strong interactions one gets 16 operators, so that the 
Lagrangian is written as 

Cf^ = -FiN^NN'^NN^NN^N ~ F2N'fT''NNW''NN^NN'fN 

-F3N^t''NN^t''NN'<t''NN''t''N - FiN'<a'NN'<a'NN''NN''N 
-F5N^aWNN''aWNN^NN''N - FeN^<7WNN''a'NN''T''NN^ N 
-F7N'<a'NN'<a'NN'<T''NN'<T''N - F^N^ a' N N'^ a' t" N N'^ N N'^ N 
-FgN'<aWNN'<aW''NN''T''NN''T''N - FwN'^ a' N N'^ a' N N'^ N'^ N 
-FiiN'^ N N'^ N N'^ N N'^ N - Fi2N'^ N N'^ N N'^ N N'^ N 
-Fi^N'^ N N'^ N N'^ N N - FuN'^ N N'^ a^T^ N N'^ N N'^ a' N 
-Fi5e^^^e''^''N'^u'T''NN'^ah^NN'^a^T''NN^N - FiQe''^t''^''N'^a'T''NN'^uh^NN'^a^NN'^T''N 

16 

^ -^^^.or- (19) 

1=1 

However, after using all possible Fierz rearrangements (cfr. Appendix VK\ . as was done in the ZN sector, it turns out 
that they are all equivalent, leaving only one independent 'iN contact operator. Using cutoff functions depending 
only on k^, e.g. choosing F(k^; A)i^(k^; A)F(k^; A), the 4-body contact potential is local in coordinate space, 

V^N^F J2 

ZQ{ri.i)Zo{rik)ZQ{rii), (20) 

where F is the single accompanying LEG. This allows to decouple, to a certain extent, the 3iV sector from the 4A^ 
sector, if one is willing to fix all the LECs from the data. For instance, one can in principle adjust the subleading 3N 
potential in the A = 3 systems without worrying much about the consequences for the a particle binding energy: the 
latter could always be reproduced by adjusting the 4A^ contact term. 
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VII. CONCLUSIONS 

In the present paper we have derived a 3N potential from the minimal form of the two-derivative 3A^ contact 
Lagrangian. This potential has 10 unknown LECs and, with a particular choice of the cutoff, can be put in a local 
form. In addition we have shown that the leading AN contact Lagrangian consists of only one operartor. 

It should be stressed that these terms start to contribute at N4L0 and are therefore beyond the accuracy of the 
presently developed nuclear interactions. In particular, a complete EFT calculation should also address the problem 
of the NN interaction at the same order. Notice that the four-derivatives 2N operators are already part of the 
N3L0 interaction while the six-derivative ones would start to contribute at N5L0, and have not been considered 
in the literature so far. Nevertheless, despite respresenting only part of the N4L0 interaction, the terms derived 
here could play an important role in the accurate description of the three-nucleon systems, since they are completely 
unconstrained by symmetries. This is at variance with the 3A^ N3L0 interaction, which contains no free parameters 
[l^ll^]. Moreover, the same terms will also appear in the pionless version of the effective theory as next-to-leading 
37V interaction and leading AN interaction. 

The utility of the derived potential could be the following: After a sensitivity study of the 10 subleading terms of 
the 3N potentials, the corresponding LECs can be determined from a fit of several polarization observables in N-d 
scattering at low energies. There are well established discrepancies between theoretical predictions and experimental 
data in some of the polarization observables as for example the vector analyzing powers Ay and iTn and the vector 
analyzing power T21 in elastic scatterin g [2 8.-^1, and in several unpolarized breakup cross sections (30l. |3]|. Our 
expectation, based on the results of Ref. [25[, is that some of the 10 operators of the 3A^ potential will have sufficient 
sensitivity to fix these discrepancies with reasonable values of their LECs. At the same time the 4N contact term 
will be used to reproduce the *He binding energy. Finally, the predictions of the derived potentials will be tested 
in the description of the four-body scattering states. In fact, also in this case there exist large and still unexplained 
discrepancies between theory and experiment [33 - l35| . Studies along these lines are in progress. 

Appendix A: Linear relations from Fierz's constraints 

In this Appendix we list the linearly independent relations among the 146 operators of Table HI The following ones 
are obtained by applying Fierz's reshuffling of indeces of particles 1-2: 



2o25 


=076 + 2077 - 3078 + 2o8i - 2o82 + Os4 ^ Ose 




(Al) 


026 


=080 - 082 + 2o85 - 2086 




(A2) 


2032 


=3080 - 3082 + 6o85 - 6086 - O120 + O122 - 2oi25 + 2oi26 




(A3) 


8035 


= - 2O3 - 20i2 + O34 + O35 + O54 + 056 + 0128 + O129 + Oi32 -|- O133 




(A4) 


2057 


= - 076 + 078 - 084 + 086 




(A5) 


058 


=082 — O80 




(A6) 


4062 


= 2oi5 - 2oi6 + 2oi9 - 2020 - O37 -1- 038 - O40 + O41 -1- 046 - O47 + O50 


- O51 - 2oi38 






+ 2oi39 + 2oi42 - 2oi43 




(A7) 


2063 


= — O80 + 082 — O116 -I- O118 + O120 — O122 — O124 + Oi26 




(A8) 


2064 


= — 3o80 + 3082 + O120 — O122 




(A9) 


4065 


=034 — O35 + O54 — O56 + Oi28 — O129 — O132 + O133 




(AlO) 


4068 


=3o34 - 3035 - O54 + 056 + 3oi28 - 3oi29 + O132 - O133 




(All) 


2083 


=079 + 081 + 083 + 084 




(A12) 


4087 


=042 + O44 — O53 — 056 + Oi30 + O133 — O134 — O135 




(A13) 


4088 


=3042 - 044 - 3053 + 056 + 3oi3o - O133 - 3oi34 -|- O135 




(A14) 


4oi23 


=3079 -1- 308l + 3083 + 3o84 - Ong - O121 - 0x23 - O124 




(A15) 


4oi24 


=9o79 - 308l + 9o83 - 3o84 - 3oiig -t- O121 - 3oi23 + O124 




(A16) 


80125 


=3026 - O32 -1- 3058 - 064 + 6O80 + 6082 + 6035 -|- 6036 - 2oi20 - 2oi22 


— 2oi25 — 2oi26, 


(A17) 



while Fierzing on particles 1-3 we get 



4oi = - Oi - O3 - 021 - 023 



(A18) 
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402 = - 02 - O3 - 022 - 024 - 065 " O74 (A19) 

403 = - 3oi + 03 - 3o2i + 023 (A20) 

404 = - 04 - 06 - 013 - O16 - O91 - O94 (A21) 

405 = - O5 - 06 - 014 - 015 + 028 " 066 " 069 - O92 - O93 (A22) 

407 = - O7 - O9 - 017 - O20 - Oi03 - O106 (A23) 

408 = - 08 - O9 - O18 - 019 - 029 - 067 " O70 - Oi04 - Oi05 (A24) 
40io = - Oio - O12 - O21 - 024 - 087 - O90 (A25) 
4oi3 = - Oi — 03 - 013 - 015 - 017 - 019 + O21 + O23 - O79 - 082 + 083 + 086 (A26) 

2025 =014 — 016 — 018 + O20 + O80 — 081 + 084 — 085 (A27) 

2026 = - 014 + O16 + O18 - O20 + O80 - 081 + 084 - 085 (A28) 
027 =011 - O12 - O22 + 023 (A29) 

2028 = - 08 + O9 + Oil - O12 + O18 - 019 - O22 + O23 + 076 " O77 - O108 + Oi09 - O112 

+ 0113 - 0116 + 0117 (A30) 

029 = - 08 + 09 + 018 - 019 (A31) 

2030 = - 08 + O9 + On - O12 + O18 - O19 - O22 + O23 - 076 + O77 + O108 - O109 + O112 

- 0113 + 0116 - 0117 (A32) 

2031 =05 - 06 - 08 + O9 - 014 + Oi5 + O18 - 019 + O76 - O77 - O116 + O117 + O120 - O121 

+ 0124 - 0125 (A33) 

4033 = — 033 — 035 — 053 — 055 (A34) 

4034 = — 034 — 035 — 054 — 056 + 065 + O71 (A35) 

4035 = - 3033 + O35 - 3053 + O55 (A36) 

4036 = - 036 - 038 " O45 - O48 - Oio7 - O109 (A37) 

4037 = - O37 - 038 - 046 " O47 + Oeo + 066 + O73 - O108 - Ono (A38) 

4039 = - 039 - 041 - 049 - O52 - 095 - 097 (A39) 

4040 = - O40 - O41 - O50 - O51 - 061 + 067 + O72 - 096 " O98 (A40) 
4042 = - O42 - O44 - O53 - 056 + 087 + 089 (A41) 
4045 = - O33 - O35 - O45 - O47 - O49 - O51 + O53 + O55 + O75 + O77 - 083 - 085 (A42) 
2O57 =046 - 048 - O50 + O52 - 076 + O78 - 084 + 086 (A43) 

059 =043 — 044 — 054 + 055 (A44) 

2O60 = — 040 + 041 + 043 — 044 + 050 — 051 — 054 + 055 — O80 + 082 — O92 + O94 + OlQO 

- 0102 + 0120 - 0122 (A45) 
061 = - 040 + 041 + 050 - 051 (A46) 

2062 = - 040 + O41 + O43 - 044 + O50 - O51 - O54 + O55 + O80 - 082 + O92 - O94 - Oioo 

+ 0102 - 0120 + 0122 (A47) 

2065 =034 - O35 + O54 - 056 (A48) 

2066 =037 - 038 + 046 - 047 + O108 - Ono (A49) 

2067 =040 - O41 + O50 - O51 + Ogg - Ogs (A50) 

2068 =043 - O44 + O54 - O55 - 088 + Ogo (A51) 

2069 =034 - O35 + 046 - 048 + O50 - O52 - O54 + O56 - O76 + O78 + 034 - 086 (A52) 

2070 =034 - O35 + 046 - 048 + O50 - O52 - O54 + 056 + 076 " O78 - 084 + 086 (A53) 

2071 =3034 - 3O35 - O54 + 056 (A54) 

2072 =040 - O41 + O50 - O51 - 096 + 098 (A55) 

2073 =037 - 038 + 046 - 047 - 0108 + Ono (A56) 

2074 =043 - 044 + 054 - 055 + 088 - 090 (A57) 

4075 =045 + O47 - O49 - O51 - O75 - O77 - 083 - 085 (A58) 

4076 =046 + 048 - O50 - O52 - O57 - O58 - 069 + O70 - O76 - O78 - 084 - 086 (A59) 
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4077 =3045 - O47 - 3049 + O51 - 8075 + O77 - Sogs + OgS (A60) 

2079 = - 079 - 081 (A61) 

2080 = - O57 - O80 - 082 (A62) 
4083 = - O45 - O47 + O49 + O51 - O75 - O77 - 083 - 085 (A63) 

2087 =042 + O44 - O53 - 056 (A64) 

2088 =043 + O44 - O54 - 055 - 068 + O74 (A65) 

2089 =3042 - O44 - 3053 + 056 (A66) 
409i = - 039 - 041 + 042 + 044 + 049 + 052 - 053 - 056 - 079 — 081 — 091 — 093 + 099 

+ Oloi + 0119 + 0121 (A67) 

4093 = — 8039 + 041 + 3042 — 044 + 3049 — O52 — 3053 + O56 — 3079 + 081 — 3091 + O93 

+ 3o99 - oioi + 3oii9 - 0121 (A68) 

2095 = - 039 - 041 + 049 + 052 (A69) 

2096 = - O40 - O41 + O50 + O51 + 067 - O72 (A70) 

2097 = - 3039 + 041 + 3049 - 052 (A71) 

4O99 = - O39 - O41 + O42 + O44 + O49 + O52 - O53 - 056 + O79 + 081 + O91 + O93 - O99 

- oioi - 0119 - 0121 (A72) 

40i03 = - 036 - 038 + O42 + O44 + O45 + O48 - O53 - O56 + O79 + 081 - O103 - Oi05 + Om 

+ 0113 - 0119 - 0121 (A73) 

4oi05 = — 3036 + 038 + 3042 — O44 + 3045 — O48 — 3053 + O56 + 3079 — 081 — 3oi03 + Ol05 

+ 3oiii - 0113 - 3oii9 + 0121 (A74) 

20107 = - 036 - 038 + O45 + 048 (A75) 

20108 = - O37 - O38 + 046 + O47 + 066 - O73 (A76) 

20109 = - 3036 + 038 + 3o45 - O48 (A77) 

40119 = - 3079 - 3o8i + 0119 + 0121 (A78) 

40120 = - 3O57 + 063 - 3O80 - 3082 + O120 + O122 (A79) 

40121 = - 9o79 + 3o8i + 3oii9 - 0121 (A80) 

40122 =3o57 - 063 - 3080 - 3082 + O120 + O122 (A81) 
I6O137 = - 2oi3 - 2oi4 - 2oi7 - 2oi8 + 036 + 038 + O39 + O41 + 8045 + 4o46 + O48 + 5o49 

+ 4050 + O52 + 2095 + 2097 + 4099 + 4oioo + 2oi07 + 2oi09 + 4oiii + 4oii2 

- 20i37 - 20i40 - 20i41 - 20i42 (A82) 
I6O140 = - 6013 + 20i4 - 6O17 + 20i8 + 3036 - O38 + 3O39 - O41 + I5O45 - 4046 - 048 

+ 15049 - 4050 - O52 + 6095 - 2097 + 120g9 - 4oioo + 60107 " 2oi09 + 12oiii 

- 4oii2 - 60137 + 2oi4o - 60141 + 2oi42, (A83) 

and on particles 2-3, 

025 =076 - 077 (A84) 

028 =05 - 06 - Oi4 + 015 (A85) 
2029 =05 — 06 — Oil + 012 — 014 + O15 + 022 — 023 + O76 — O77 + O96 — O97 + OlQO — OlQl 

- 0116 + 0117 (A86) 

2O31 =3076 - 3077 - O116 + 0117 (A87) 

O57 =082 - O80 (A88) 

O60 =037 - 038 - 046 + O47 (A89) 

2061 =037 - 038 - O43 + O44 - 046 + O47 + O54 - O55 - O80 + 082 + O104 - O106 - O112 

+ 0114 + 0120 - 0122 (A90) 

2063 = - 3O80 + 3082 + O120 - O122 (A91) 

2065 = - O2 + O3 - O22 + 024 (A92) 

2066 = - O5 + 06 - 014 + 015 - O92 + O93 (A93) 
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2067 = - Os + Og - Ois + Oi9 - Oio4 + Oi05 (A94) 

2068 = - On + Oi2 - 022 + 023 - 088 + 089 (A95) 

2069 = - 05 + 06 - 014 + 015 + O92 - O93 (A96) 

2070 = - OS + 09 - 018 + 019 + oio4 - 0105 (A97) 
2on = - On + 012 - 022 + 023 + 088 - 089 (A98) 

2072 = - 02 + 03 - Oi4 + O16 - O18 + O20 + O22 - 024 + O80 - 081 - 084 + 085 (A99) 

2073 = - 02 + 03 - Oi4 + O16 - O18 + O20 + O22 - 024 - O80 + 081 + 084 - 085 (AlOO) 

2074 = - 302 + 303 + 022 - 024 (AlOl) 

2075 = - 075 - 078 (A102) 

2076 = 025 - 076 - 077 (A103) 
4o79 = - O13 - O15 + O17 + O19 - 079 - 082 - 083 - 086 (A104) 
4082 = - 3oi3 + O15 + 3oi7 - 019 - 3o79 + 082 - 3o83 + 086 (A105) 
4084 =0i4 + O16 - O18 - O20 + O25 + 026 - O72 + O73 - O80 - 081 - O84 - O85 (AIO6) 

2o87 = - oio - 012 + 021 + 024 (A107) 

2090 = - 3oio + 012 + 3o2i - 024 (A108) 

2091 = — 04 - 06 + oi3 + 016 (A109) 

2092 = - 05 - 06 + 014 + 015 - 066 + 069 (AllO) 

2094 — — 3o4 + 06 + 3oi3 - 016 (Am) 

4095 = — 04 - 06 + Olo + 012 + 013 + 016 — O21 — O24 - O75 - O78 - O95 - O98 - O99 

- 0102 + 0115 + 0118 (A112) 
4o98 = - 3o4 + 06 + 3oio - 012 + 3oi3 - 016 - 3o2i + 024 - 3075 + 078 - 3095 + 098 

- 3099 + 0102 + 3oii5 - 0118 (A113) 

20103 = - 07 - 09 + oi7 + 020 (A114) 

20104 = - Os - O9 + O18 + Oi9 - 067 + O70 (A115) 

2oio6 = - 3o7 + 09 + 3oi7 - 020 (A116) 

40115 = - 3075 - 3078 + O115 + O118 (A117) 

40116 =3025 - O31 - 3076 - 3o77 + O116 + O117 (AII8) 

40117 = - 3025 + O31 ~ 3076 - 3077 + O116 + O117 (A119) 

40118 = - 9o75 + 3o78 + 3oii5 - 0118 (A120) 
4oi27 = - 0127 - 0129 - 0134 - 0136 (A121) 
2oi28 = - 0128 - 0135 {A122) 
4oi29 = - 3oi27 + 0129 - 3oi34 + O136 (A123) 
2oi30 = - O130 - 0132 (A124) 
4oi3i = - 2oi3i - 2oi33 + 0144 (A125) 
2oi37 = - O137 - 0139 (A126) 
4oi38 = - 2oi38 - 2oi4o + 0145 (A127) 
4oi4i = - 0127 - 0129 + 0134 + O136 - 2oi4i - 2oi43 (A128) 
2oi42 = - 0128 + 0135 - 20i42 (A129) 

0144 =2oi3i - 2oi33 (A130) 

0145 =2oi38 - 2oi4o (A131) 

0146 = - 0131 + 0133 + 0138 " O140. (A132) 

The above relations have been selected according to a criterium of simplicity. Indeed, quite lengthy combinations may 
arise: for instance, exchanging particles 1-2 in the operator 093 gives rise to a combinations of 63 different operators 

of Table ID All other relations are linearly dependent on the selected ones. The algebraic manipulations have been 
authomatized using the program FORM |36| . 

We now list 15 independent linear relations among the 16 operators defined in Eq. ((T9)l . out of 96 linear relations 



from Fierz's reorderings of 6 pairs of the four particles: 
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All 16 operators in Eq. ^T9\\ are therefore proportional to each other. 
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